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We study vortical states in a Bose-Einstein condensate (BEC) filling a cigar-shaped trap. An ef- 
fective one-dimensional (ID) nonpolynomial Schrodinger equation (NPSE) is derived in this setting, 
for the models with both repulsive and attractive inter-atomic interactions. Analytical formulas for 
the density profiles are obtained from the NPSE in the case of self-repulsion within the Thomas- 
Fermi approximation, and in the case of the self-attraction as exact solutions (bright solitons). A 
crucially important ingredient of the analysis is the comparison of these predictions with direct 
numerical solutions for the vortex states in the underlying 3D Gross-Pitaevskii equation (GPE). 
The comparison demonstrates that the NPSE provides for a very accurate approximation, in all the 
cases, including the prediction of the stability of the bright solitons and collapse threshold for them. 
In addition to the straight cigar-shaped trap, we also consider a torus-shaped configuration. In 
that case, we find a threshold for the transition from the axially uniform state, with the transverse 
intrinsic vorticity, to a symmetry-breaking pattern, due to the instability in the self-attractive BEC 
filling the circular trap. 

PACS numbers: 03.75.Lm,03.75.Kk,03.75.Hh 
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I. INTRODUCTION 

An important setting for experimental and theoretical 
studies of dynamical phenomena in Bose-Einstein con- 
densates (BECs) is provided by nearly one-dimensional 
(ID) "cigar-shaped" traps, which assumes tight confine- 
ment in the transverse plane, allowing to unravel the dy- 
namics along the longitudinal axis. The use of this geom- 
etry helped to achieve famous experimental results, such 
as the creation of single [l| and multiple ^ bright soli- 
tons in the condensate of several thousand ^Li atoms. 
In those experiments, the strength of the interaction be- 
tween atoms was controlled and made weakly attractive 
by means of the Feshbach-resonance technique. In the 
condensate of ^^Rb atoms trapped in a similar geometry, 
nearly 3D solitons were found in a post-collapse state [1] . 

It is natural to derive effectively one-dimensional (ID) 
equation(s) for the description of this experimentally rel- 
evant setting, starting from the full 3D Gross-Pitaevskii 
equation (GPE). The reduction of the 3D equation to a 
ID form was performed, under various assumptions, by 
means of sundry methods [J|-H- ^^ ^^^ simplest approx- 
imation, the difference of the eventual equation from a 
formal ID variant of the GPE (with the underlying cu- 
bic nonlinearity) is represented by an additional quintic 
term, whose sign is always self-attractive, irrespective of 
the sign of the cubic nonlinearity [6| . A more consistent 
approach to the derivation of the ID equation postulates 
the factorization of the 3D mean-field wave function into 
the product of the ground state of the 2D harmonic os- 
cillator in the transverse plane, and slowly varying axial 
(ID) wave function. Using the variational representation 
of the underlying 3D GPE, this approach ends up with 



the nonpolynomial Schrodinger equation (NPSE) for the 
axial wave function [1, & [lOl, [HI- The NPSE has been 
extended to investigate the Tonks- Girardeau regime [12|. 
the two-component BEC pJl , and also transverse spatial 
modulations |14| . 

Another physically interesting modification of the 
above setting is that when the BEC trapped in the cigar- 
shaped geometry is lent vorticity (assuming that the cor- 
responding vector of the angular momentum is directed 
parallel to axis of the cigar-shaped trap), as proposed in 
Ref. [l5|. In the experiment, the vorticity may be im- 
parted on the condensate by a helical laser beam shining 
along the axis of the trap. A natural problem for the the- 
oretical analysis, which is considered in the present work, 
is the derivation of a modification of the effective one- 
dimensional NPSE that takes into regard the intrinsic 
vorticity. As we demonstrate in Section II, such a modi- 
fied ID equation is not drastically different from its coun- 
terpart derived before [^, Q for the zero- vorticity states. 
However, a really nontrivial issue is to compare basic dy- 
namical states predicted by this effective equation with 
direct numerical solutions of the underlying GPE in three 
dimensions. In this work, we perform the comparison 
separately for the self-repulsive and self-attractive BEC 
(in Sections III and IV, respectively). In the former case, 
the relevant states are of the Thomas-Fermi (TF) type, 
while in the latter model (self-attraction) these are bright 
solitons. We demonstrate that the modified NPSE pro- 
vides for very good approximation for all these states, 
including the prediction of their stability (and, in par- 
ticular, of the collapse threshold for solitons in the self- 
attraction model). In addition, in Section V we consider 
a different geometry, when the long quasi-lD trap is bent 



and closed into a ring (torus), maintaining the intrinsic 
vorticity in the transverse plane. In that case, we find 
that the modified NPSE accurately predicts the shape 
and stability of a possible soliton, as well as the dclo- 
calization transition to the state of the BEC uniformly 
filling the toroidal trap. 



II. THE NONPOLYNOMIAL SCHRODINGER 

EQUATION (NPSE) FOR MATTER- WAVE 

VORTICES 

The description of a dilute BEC, confined in the axial 
direction by a generic potential V{z), and in the trans- 
verse plane by the harmonic potential with frequency lu_i_, 
is based on the fundamental GPE in three dimensions 
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where V'(r, t) is the macroscopic wave function of the con- 
densate normalized to unity, and g = 2Na/aj_, with N 
the number of atoms, a± ~ ^Jfi/ {mu}±) the length of 
the transverse harmonic confinement, and a the scatter- 
ing length of atomic collisions. In Eq. (IT|), the length 
and time units are a^ and wj , and the energy unit is 
fiLj^. This equation can be derived from the Lagrangian 
density: 
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In this work, we study the existence and properties of 
matter- wave states with vorticity in the transverse plane, 
which correspond to solutions of the form 



ip{r,t) = <i>(r, z,t) e 
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where S is the integer vorticity quantum number, r and 
6 being the polar coordinates in the [x, y) plane. On the 
substitution of expression ([3]), Eq. ^ takes the form: 
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We solved Eq. ^ numerically by using the finite- 
difference Crank-Nicholson predictor-corrector method 
with the cylindric symmetry described in Ref. [16:] . The 
Crank-Nicholson algorithm with imaginary time is used 
to find the stationary wave function of the GPE with a 
fixed vorticity S. The algorithm provides for fast con- 
vergence to a finite-size wave function with finite energy. 
In the case of instability due to collapse {g < 0) the al- 
gorithm quickly approaches a delta peak with (negative) 
diverging energy. Note that, in the free space, vortices 
are usually destroyed by the instability against azimuthal 



perturbations Il7l |, w hich split them into a set of sepa- 
rating fragments [l8|, but in our setting this instability 
is suppressed by the tight transverse confinement. 

To gain some analytical insight and to reduce the prob- 
lem to 1-D, we apply a variational approach. In the case 
of a cigar-shaped geometry, it is natural to extend the 
variational representation 0, @] of the 3D GPE which led 
to the NPSE for the axial wave function in the case of 
no vorticity to the present situation. Thus we adopt the 
following ansatz for the vortex state described by Eq.Q: 
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where a{z, t) and /(z, t) which account for the transverse 
width and the amplitude of the vortex are to be deter- 
mined variationally. By inserting this ansatz into the 
Lagrangian density ^ and performing the integration 
over X and y, we derive the effective Lagrangian density. 
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where the effective nonlinearity strength for the vortex 
state is 
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Note that the effective Lagrangian density ^ does not 
contain time derivatives of a{z,t). Expression ([6]) gives 
rise to the corresponding Euler-Lagrange equations. 
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where a' = da/dz. Neglecting this derivative, Eqs. ([7]) 
and dHl) can be combined into: 
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which is a new version of the nonpolynomial Schrodinger 
equation (NPSE) for the vortical state tightly confined 
in the transverse plane. In the case of 5 = 0, Eq. ([9|) 
reduces to the ordinary NPSE, that has found various 
applications [a. [l4|. 

The eigenfunctions of the equation 
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provide the stationary solutions of Eq. ([9]) of the form 
f{z,t) = (t){z) e-^^*. Obviously, Eqs. (© and 1^ can 
be transformed into their counterparts for 5 = by 
substitutions g — > gs, /x -^ /^/(-S + 1), t ^ (5*+ l)i, 
z -^ VSTTz, and V ^ V/{S + 1). Therefore, the re- 
sults obtained in earlier works in the framework of the 
NPSE with S — can be easily apphed to Eqs. ([9]) and 
(fTU)) . A really nontrivial issue is to compare these results 
with those produced by numerical integration of the full 
three-dimensional GPE for vortical states. 



III. REPULSIVE NONLINEARITY 

In the case of a repulsive inter-atomic interactions (i.e. 
g > 0) one may readly generalized the TF approximation 
developed in Ref. 10], for 5 = to the case with vor- 
ticity iS*. Thus, neglecting the kinetic-energy term, i.e., 
the second derivative, in Eqs. (jTU]) one gets the analyti- 
cal expression for the stationary axial probability density 
p(z) = |(/'(z)P in the form: 
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where the effective local chemical potential is 
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Notice that generally expression (fTTjl does not represent 
a soliton. Equations (|11|) and (fT2|) are a straightforward 
gen eralization of the TF approximation developed in Ref. 
|l3 for the repulsive BEG with S = 0. 

To test the accuracy provided by the NPSE in this 
case, in Fig. [1] we plot the axial and radial probability 
densities of the repulsive BEG with vorticity S, defined, 
respectively, as 
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and obtained from the numerical solution of the station- 
ary version of Eq. ^ . They are compared with predic- 
tions of the stationary NPSE (jlOp , as well as with an- 
alytical result pT|) produced by the TF approximation. 
We choose g ~ 20, and take V{z) — z^/2. This choice 
of the axial potential implies that the full potential is 
isotropic. Thus, we are testing the NPSE in a geometry 
far from the cigar-shaped one where it is expected to be 
very accurate, as it was verified in the case 5 = 0. 

Figure [1] shows that, for all S, the agreement between 
the 3D GPE and the NPSE (sohd and dashed hues, 
respectively) is extremely good in the axial direction. 
There are some differences in the radial probability den- 
sity, p{r). Remarkably, the TF version of the NPSE, i.e., 
Eq. (fTTj) (which corresponds to the dotted-dashed line 
in the figure) also provides for very good accuracy of the 
axial probability-density distribution in the core of the 
confined space. 
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FIG. 1: (color online). Axial p{z) and radial p{r) probabil- 
ity densities of the state with vorticity S in the waveguide, 
with strength of self-repulsion g = 20, and the harmonic 



potential applied also in the axial direction, V{z) 



72. 



Solid lines: numerical solution of the stationary 3D Gross- 
Pitaevskii equation (GPE); dashed lines: numerical solution 
of the effective one-dimensional NPSE, Eq. (|10|l : dotted lines: 
obtained the analytical Thomas-Fermi (TF) approximation 
obtained from the NPSE, Eq. (fTTj) . Dot-dashed lines: the 
TF approximation for the 3D GPE. 



The comparison between the full GPE and the effec- 
tive NPSE is further presented in Table 1, which displays 
values of the ground-state energy of the repulsive BEG 
with vorticity S, as obtained from the three-dimensional 
GPE (Egpe), and from the NPSE (Enpse)- The former 
energy is calculated as 

(13) 
where $(r) is the ground-state stationary solution of the 
GPE in the 3D setting. The energy Empse is also ob- 
tained from Eq. ([13]), but with $(r) corresponding to 
the ground-state stationary solution of the NPSE. The 
third column in the table shows that the relative error 
AEnpse/Egpe is quite small, and moreover that de- 
creases with the increase of S. 
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AEnpse/Egpe 


AEtf/Egpe 
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3.0729 
3.6704 

4.4878 


0.0194 
0.0071 
0.0070 


0.1855 
0.2558 
0.2295 



Table 1. Repulsive Bose-Einstein condensate with nonlin- 
earity strength g — 20 and vorticity S, in axial trapping 
potential V{z) — z'^/2. Egpe represents the ground-state 
energy as found from the full 3D GPE. AEnpse/Eqpe and 
AEtf/Egpe are the relative differences of Egpe with re- 
spect to the ground state energy of NPSE and TFGPE wave 
functions. 

In addition, we have also calculated the size of the 
term {S + l)(T'^/cr^ in the effective Lagrangian ©, that 



was neglected in the derivation of the NPSE. The result is 
that the neglected term is very small indeed, being always 
less than 0.5% of the energy produced by the NPSE. 

For the sake of completeness, we have also calculated 
the axial probability density obtained by the application 
of the TF approximation to the full 3D GPE (TFGPE) 
and subsequent integration over x and y. It is given by 
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where /i(z) = /i — V{z). The left panels of Fig. [l]show 
that the results provided by this TFGPE approximation 
are slightly better than those obtained by the same ap- 
proximation applied to the NPSE, but the full solution 
of NPSE is better than the TFGPE approximation, as 
shown in the last column of Table 1 and slso in the left 
panels of Fig. 1, where the NPSE curves practically coin- 
cide with their counterparts produced by the numerical 
solutions of the GPE in 3D. 



IV. ATTRACTIVE NONLINEARITY 

Proceeding to the attractive inter-atomic interactions, 
i.e., negative g (and gs), it is relevant to recall that, for 
V{z) — and 5 = 0, a family of bright-soliton solutions 
to Eq. ^ was constructed in Refs. ,5, _9J. In the case of 
5 7^ 0, one is dealing with bright vortex solitons. In a nu- 
merical form, they were studied in Ref . [l9| . A Gaussian- 
based variational approach was developed in Ref. [15| . 
and, in a brief form [which resulted in an effective poten- 
tial for the vortex soliton in the case of inhomogeneous 
transverse confinement, uj± = uj±{z)] also in Ref. [8|. 
The stability of the vortex solitons with both large and 
small aspect ratios (i.e., strongly elongated ones, as we 
consider here, or pancake-shaped solitons, tightly con- 
fined in the axial direction) was studied by means of ac- 
curate numerical methods, based on the linearization of 
the GPE with respect to small perturbations, in Refs. 

m- 

Here we use the NPSE in the form of Eq. (fTU)l and 
compare its results with those from the full 3D GPE. To 
simplify the notation, we set 
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By imposing vanishing boundary conditions at infinity, 
(l){z) ^ as z — > ±00, we find soliton solutions in an 
implicit analytical form. 
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The soliton family is then characterized by the depen- 
dence of 7s on ^s- Eq- (0, implies that solitons do 
not exist if the nonlinearity strength exceeds the criti- 
cal value, 7^"^ = 4/3, corresponding to /^s = 1/2. For 
stronger nonlinearities the NPSE ^ predicts a longitu- 
dinal collapse, that is a quantum tunneling to a high- 
density state |2(?| where atoms quickly evaporate due to 
three-body recombination (3|. In terms of the physical 
parameters, a < 0, a±^ and N (number of atoms), we 
conclude that the collapse of the bright vortex soliton 
takes place in region 
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In the opposite limit of weak nonlinearity, 75 
soliton solution (fT4|) takes the ordinary form. 



(16) 



0, exact 



Hz) = y^sech fl^sVS + Tz) . 

As a final remark, we note that the inversion of Eq. 
()15p provides two branches of solutions for /is (7s)- In 
numerical simulations of Eq. ([5]), only the one satisfying 
the condition d^s/d'ys < turns out to be stable, in 
precise agreement with the pre diction of the Vakhitov- 
Kolokolov stability criterion [2l| . 

In Fig. [2] we plot the axial probability density p{z) of 
the soliton with vorticity S*, with no axial external poten- 
tial {V{z) = 0), and for 5 = —1. The figure shows that 
the radial probability density p{r) produced by the NPSE 
is virtually indistinguishable from that generated by the 
3D GPE. This can be understood since the system with 
g — —\ becomes effectively quasi-one-dimensional. The 
figure also shows that, for all S, the agreement between 
the 3D GPE and the NPSE (solid and dashed fines) is 
excellent also in the axial direction. Thus, we conclude 
that the NPSE predicts the bright vortex solitons with a 
very high accuracy. 



V. VORTEX BRIGHT SOLITONS IN A RING 

Repulsive BECs in a ring-shaped (toroidal) trap are 
now available to the experiment [24! . Self-attractive BEG 
in a ring has not been experimentally created yet, but it 
may also be a promising setting, interesting from a phys- 
ical point of view, since a quantum phase transition from 
a uniform state to a bright soliton has been predicted in 
Refs. [2J] . This prediction is based on the mean-field and 
beyond-mean-field numerical analysis for the ID Bose gas 
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FIG. 2: (color online). Axial p{z) and radial p{r) probability 
densities in the soliton state with vorticity 5", in the case of 
the self-attraction with strength g = —1 and V{z) = (no 
axial potential). Solid lines: obtained from the full 3D GPE. 
Dashed lines: produced by the NPSE. 



with contact interactions and periodic boundary condi- 
tions. Recently, self-attractive BECs in a 3D ring-shaped 
trap were considered, taking into account the transverse 
structure of the trapped condensate. In this way, new fea- 
tures in the model's phase diagram have been predicted, 
both at zero [2J| and finite [25| temperatures. 

The analysis presented in Refs. [21, [2^ can be ex- 
tended to the bright sohtons with intrinsic vorticity by 
using the above-mentioned scaling that makes Eq. (fTU]) 
equivalent to its counterpart with S — 0. As shown in 
Ref. 24] , one can model the BEC in the tight 3D ring by 
using Eq. (|10p with natural periodic boundary conditions 
(b.c), (j){z + L) = (j){z), with -L/2 < z < L/2, where 
L — 2ttR is the length of the ring with radius R (here, 
we do not assume any global vorticity imposed along the 
closed ring). 

Equation (fTO]) subject to the periodic b.c. always ad- 
mits the axially uniform solution, cj) = 1/^/L. Further, 
in Ref. [23 it was shown that this uniform state, with 
5" = 0, is energetically and dynamically stable in the at- 
tractive BEC only for a finite range of parameters. In 
the case of S* 7^ 0, the stability condition is modified, by 
the above-mentioned scaling transformation, to 
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In terms of physical parameters, condition (|17p takes the 
following form, cf. condition (|16|) for the onset of the 
collapse in the rectilinear trap: 
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FIG. 3: Axial probability density p{z) in the state with trans- 
verse vorticity S, in the ring-shaped waveguide, i.e., with the 
periodic boundary conditions in z. Self-attraction strength 
is <; = —1, and the length of the ring is 40. The results are 
obtained from the NPSE. 



Note that in this inequality both the scattering length a 
and L are expressed in dimensional units. 

In Fig. [3|we plot several profiles of the axial probability 
density p{z), as obtained from the NPSE, Eq. ((9|), with 
vorticity S and g = —1 [with no external axial potential: 
V{z) = 0]. Periodic boundary conditions are imposed, 
with period (length of the ring) L — 40. We have checked 
that, as well as in Fig. [2l the agreement between results 
produced by the 3D GPE and NPSE is very good for all 
S (we display only the profiles generated by the NPSE 
since they are indistinguishable in the figure from their 
counterparts obtained from the 3D equation). Note that, 
contrary to the case shown in Fig. ^ in the present case 
the axial probability density is uniform for S = 2, due to 
the delocalization transition, the threshold for which is 
very accurately predicted by Eq. (fT8| . 



VI. CONCLUSIONS 

In this work, we have presented a systematic analysis 
of self-repulsive and self-attractive matter- wave patterns 
with the intrinsic vorticity which fill a straight or toroidal 
cigar-shaped trap. The ID NPSE was derived from the 
underlying 3D GPE, as a modification of the previously 
derived ID equation for the states without vorticity. In 
the cases of the self-repulsion and self-attraction, the 
modified NPSE predicts, respectively, the axial density 
profiles for the TF states and bright solitons, both in 



an analytical form. The comparison of these predictions 
with direct numerical solutions of the underlying 3D GPE 
for the vortex states demonstrates a very high accuracy of 
the description provided by the NPSE, which includes the 
prediction of the stability of the bright solitons, as well as 
the threshold for the onset of collapse. For the toroidal 
configuration of the trap, a threshold for the transition 
from the axially uniform vortex state to a quasi-soliton 
symmetry-breaking pattern was found too. 

This work may be naturally extended in several direc- 
tions. In particular, a challenging problem is collision of 
solitons with opposite signs of the intrinsic vorticity, +S 
and —S. This setting cannot be described by the NPSE 



(which assumes a single value of S throughout the elon- 
gated trap), hence it should be studied in the framework 
of the full 3D GPE. Another interesting generalization 
may be formation of gap solitons with intrinsic vorticity, 
in the case when the nonlincarity is sclf-repulsivc, and the 
cigar-shaped trap is combined with a periodic axial po- 
tential, like in the well-known experimental setting used 
for the creation of ordinary gap solitons [2g . 
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by GNFM-INdAM and thanks Alberto Cetoh for useful 
discussions. 
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